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ABSTRACT 

It is natural to associate optimal estimation and opti- 
mal statistical testing. In this paper continuous function 
estimation of the cumulative distribution is used to define 
two test statistics that compete with the Kolmogorov 
statistic. The first statistic, C N , is attributed to Pyke 
and the second, , is obtained by polygonalizing the sample 
distribution function. It is known that both are asymp- 
totically equivalent to the Kolmogorov statistic. Using 
the methods of J. Durbin, the small sample distributions 
are tabled as well as the critical points for significance 
levels of .20, .10, .05, .025, and .01. It is shown that 
R^ is stochastically smaller than D N , and it appears that 



C XT is also smaller than D. 
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I. 



INTRODUCTION 



The popular method for testing the simple hypothesis 

H : F = F , where F is the unknown cumulative distribution 
o o 

of the population sampled and F q is the proposed continuous 
distribution of the sample, is the Kolmogorov Djq statistic 
[Ref. 4]. The asymptotic distribution of was developed 
by Kolmogorov [Ref. 2] and tabled by Smirnov [Ref. 6]; the 
small sample distribution was tabled by Massey [Refs. 3 and 
4] and has been treated recently (among others) by Durbin 
[Ref. 1] . A modification of the Kolmogorov test has been 
proposed by Pyke [Ref. 5] and the distribution of this test, 
called C N , can be obtained from Durbin [Ref. 1]. 

The Kolmogorov D N statistic is associated with a step 
function estimation of F in the computation of the test 
statistic. This function, however, is from a vacuous class 
under the assumptions of the method, since a priorily F is 
assumed to be continuous. Of course the estimator converges 
to the continuous function F as the sample size becomes 
large without limit. 

The question of estimating the distribution F is sep- 
arate from testing H q . However, it is natural to associate 
the two and to expect optimal estimation to yield optimal 
testing methods. Two of the statistics considered in this 
paper may be viewed as being based upon continuous esti- 
mators of F. It is shown that each is asymptotically 
equivalent to D . 
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The statistics investigated are presented in section 
two. The asymptotic nature of each and the relation between 
each is shown in section three. Section four presents the 
computational method used in tabling the distribution of the 
statistics introduced in section two. Final comparisons 
are covered in section four. The tables are in the 
appendices . 
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II. BASIC DEFINITIONS, KOLMOGOROV STATISTIC, PYKE 



STATISTIC , AND A PROPOSED STATISTIC 

Let , Y 2 /. . . , Y^ be a random sample from a continuous 
population with cumulative distribution function F (y) . The 
methods presented below are for testing the simple hypothesis 

H : F = F 
o o 

The procedures discussed are based on the order statistics 
of the sample where 

-00 < — ^2 ~ * * * — < 00 



are the ordered sample points/ and the sample cumulative 
distribution is 



0 



X < x 1 



F n (x) = j/N 




< X < X 



j+1 



1 



X > X N . 



The Kolmogorov statistic for testing H q is 

d n = s “ p |f n ( x) • F o (x) 1 

and has the property of being distribution free under H q 
[R ef. 4]. This property allows for computation of the cum- 
ulative distribution of D to be done with F_ assumed to be 

N o 

uniform (0,1) without loss of generality. 
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Lemma 1 : l <_ j <jsi ^ I X j 



(j/N) | , |x. 



(j-l)/N| }. 



Proof: After transformation under H the statistic 

o 

becomes 




sup 

0<x<l 



|x 




The assertion is obvious from Figure one. 

A modification of the Kolmogorov statistic is attributed 
to Pyke [Ref. 5] where the test statistic is 

C N - l<j^N |X j - V< n +D| . 

t h 

Implicit in this method is the idea that the j order 
statistic estimates the ( ( j + l)/N) (100) percentile of F. 

Such estimators are unbiased when F q is uniform (0,1). 

The statistic C N may be viewed as being based on a 
continuous function estimator of F constructed by connect- 
ing the points (0,0); (X 1 ,1/(N+1) );...; (X^N/ (N+l) ) ; 

(1,1) with straight line segments. This is apparent since 
any point on the line segment L from (X^ , j/(N+l)) to 
(X , (j+l)/(n+l)) is a convex combination of the endpoints 

and hence the maximum separation between points on the line 
Y = X and L must occur at one of the endpoints of L. 

It is of interest to consider a sample cumulative dis- 
tribution formed by connecting the points (0,0); 

(X^ , 1/N );...; (X N ,1) with straight line segments and the 

resultant statistic 
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Z|u- 




Figure 1. Lemma 1 
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*N= S X P 



|S N (X) - F q (x) | 



where S N (x) is the sample cumulative distribution modified 
as described above. Properties of are presented in 
section three and its distribution for finite N is tabled 
in Appendix B. 
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III. ASYMPTOTIC EQUIVALENCE OF THE TEST STATISTICS 



It is well known that 

00 r _ 

Lim P[/!~ D > Z] = 2 Z (-1) exp {-2r 2 Z 2 } 

N-* 00 r=l 

[Ref. 2] and it is shown below that both/N R^ and /n C^ 
have the same limiting distribution. 

Lemma 2: For all N and all F R^ is stochastically 
smaller than , i.e., Rj^ £D N . 

Proof: By lemma one 

D N = ^ { l X j " j/N l ' l X j " (j-D/Nl) 
and it follows from the definition that 




max 

j 




j/N|} . 



The result is immediately obvious. The latter is maximized 
over a subset of those for which the former is maximized. 
Lemma 3: For all c > 0 



Lim P[/N|C N -R N |<£]=1 . 

N->°° 



Proof : 




Consider 




_ max 

j 


1*3 


- j/(N+l) | 


_ max 


|x. 


- j/N + j/N - j/(N+l) | 


j 


3 


max 

“ j 


1 x . 

1 D 


- j/N | + 1/ (N+l) = + 1/ (N+l ) 
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Thus C N - R n < 1/(N+1) 



Now consider 



R N = T l X j - 3/N 



max 



, X . - j/(N+l) + j/(N+l) - j/N 



max 
j ' D 



< -J" I X_. - j/(N+l)| + 1/ (N+ 1 ) = C N + 1/ (N+l ) 



Combining the two it is seen that 



/N |R N - C N | < /N/(N+1) 



Thus 



I*n - C N I < e] = 1 



for all N so large that 



v/N / (N+l) < e 



Lemma 4: For all e > 0 it will be shown that 



Lim P[>/N |D - Rj < e] = 1 
N-*-°° 



Proof: Let Q. T = n |X. - ( j — 1 ) /N I and note that 

N 1<0<N 1 j J ' 



Q N =l<j<N l X j - + ^ N - <3-l)/N| i ^ + 



1/N 



Thus 



v'N |Q n - Rj < 1//N . 
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Also 



°N = l<j<N { l X j " ^ /N I ' l X j " (3-l)/N|} 



max { / Q n ) 



Thus | | | max { 0 , (Q^ - R^) } | <. | - R^ | and 

it follows that 



Lim P[/N |D - Rjjl < E] = 1 



for all N so large that /N > 1/e . 

00 r-1 22 

Lemma 5: Lim P[/N R^ > Z ] = 2 Z (-1) exp {2r Z } 

ISH-oq r=l 

Proof: P[/N D n > Z] = P[/N R n + (D n - R^) > Z] 

Lim P[/NR^ + /N (D - R^) > Z] = Lim P [ /N ^ > Z] 

N-*-oo n->°o 
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IV. COMPUTATIONAL METHOD 



The method used to evaluate the cumulative distribution 
function of the statistics , C N , and is from the work 
of Durbin [Ref. 1] . A sketch of his results and the modifi- 
cations necessary for the statistics under consideration are 
presented in this section. 

"Suppose that 0 <. <. <. . . . < < 1 is an ordered 

random sample of independent observations from a uniform 
(0,1) distribution. The sample distribution function is...", 
F N (x) as defined in section two. "Let S denote the sample 
path of F^(x) as x moves from 0 to 1. In this paper we con- 
sider the probability p n (a,b,c) that S lies entirely in the 
region R between the lines ny = a+(n+c)x and ny = -b+(n+c)x, 
(a>0, b>0, a+b>0, and a-c>0)."^ 

Define : 



H 



1-6 



( 1-6 )/ 2 ! 
( 1-6 3 ) / 3 I 



P P 

1-6 -@+h 
P! 



1 

1 

1 / 2 ! 



l-@ 



P-1 



(P-1) ! 



0 

1 

1 



0 

0 

0 




0 

0 

0 

• 

1 

l-@ 



^"Durbin, J. , "The Probability That the Sample Distri- 
bution Function Lies Between Two Parallel Straight Lines." 
The Annals of Mathematical Statistics, V. 39, p. 398, 

April 1968. 
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[z] = the greatest integer in z. 

1-6 = (b-c) - [b-c] 
l-@ E (a-c) - [a-c] 

P E [b-c] + [a+c] + 1 

h E 0 if 6+@ < 1 and h E (6 + @ - 1) P if 6+@ > 1 

q e the ([b-c + l] / [b+l]) element of the matrix 
H [N+c ] 

Thus p (a,b,c) = — — ! ■ M q [Ref. 1]. 
n (N+c) 

Lemma 6: P [D^ < k/N] = p (k,k,0) [Ref.l]. 

Lemma 7: P [C N <.k/(N+l)] = p (k,k+l,l) [Ref. 1] . 
Lemma 8: P [ R^ <. k/N] = p n (k,k+l,0) 



Proof: Assume the sample cumulative distribution, 

S^(x), has stayed within the bounds shown in Figure two 
until it arrives at the point (X^,j/N). The next observa- 
tion, X j + i' must occur before Xj + ^ for non-rejection by 
the statistic. For R M to accept the sample, the next 



N 



,R 



observation must occur before Xj + ^. Since all steps are of 
height 1/N the point (X P + ^, (j+l)/N) is in the region R if 

the lower band of the equivalent statistic is shifted 
downward by 1/N. The resultant band is ny = -(b+1) + nx. 

If the point (X^ + ^, (j+l)/N) were to fall in the region W 

for either or R r . the hypothesis, H q , would be rejected. 
Hence the upper bound for the region R is unchanged. 

The computer program used to table the distribution of 
and R^ is presented in the rear of this thesis. As 
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written the parameter P is restricted to be less than or 



equal to forty, 
matrix products 



The subroutine DURMAT is used to keep the 
within the range of the IBM 360/67 used. 
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V. COMPARISON OF STATISTICS, REMARKS , AND CONCLUSIONS 



The distribution of the statistics C N and R^ are pre- 
sented for the first time in Appendix A and B. Some 
critical points for various levels of significance are 
presented in Appendix C for , C N , and R^. It is seen 
for sample sizes greater than sixty all three statistics 
have essentially the same critical points for alpha = 

.20, .10, .05, .025, and .01. 

From lemma two it is known that R^ is stochastically 
smaller than D N and Appendix C indicates that C N is also 
stochastically smaller than for the alpha values shown 
under H q . Such results are important in developing smaller 
confidence bands for distributions [Ref. 7]. 

It has been noted that the asymptotic formula of Kol- 
mogorov is correct (to three decimal places) if N>60 for all 
three statistics. The tables indicate that the limiting 
distribution is reached by earlier (i.e., N>_30) for the 
alpha values listed. It can be seen in Appendix C that for 
alpha >..05 the distribution of approaches the asymtote 
from above whereas both C N and R N approach it from below for 
all listed values of alpha. Also, and R^ are in essen- 
tial agreement for sample size greater than fifteen. 

Other statistics can be investigated empirically and 
analytically using Durbin's paper. An interesting area for 

investigation would be to relate the power functions of the 

( 

test statistics to the parameters (a,b,c). 
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APPENDIX A, 



P( C < K/ ( N+l ) ) 

N 



N = 5 1C 15 



K 



l 


.27006 


.03550 


.00417 


2 


.83719 


.46822 


.23292 


3 


.98666 


. 83796 


.63607 


4 


.99999 


. 96862 


. 87604 


5 


1.0 


.99626 


.96838 


6 




.99974 


.99403 


7 




.99998 


.99917 


a 




.99999 


.99991 


9 




. 99999 


. 99998 


10 




1.0 


.99998 


11 

12 


N 


= 40 


5 


50 


K 

3 


. 10016 


.06674 


.04420 


4 


.33714 


.26850 


.21257 


5 


.58829 


.51376 


.44601 


6 


.77448 


.71374 


. 65384 


7 


.88906 


. 8^738 


.80305 


8 


.95078 


.92581 


. 89700 


5 


.98030 


.96708 


.95073 


10 


.99288 


.98666 


.97796 


11 


.99766 


.98666 


.99097 


12 


.99929 


.99505 


.99657 


13 

14 


N 


= 75 


8C 


35 


K 

3 


.00532 


.00345 


.00224 


4 


.06238 


.04843 


.03753 


5 


.20743 


.17657 


. 15001 


6 


.39782 


.35731 


.32030 


7 


.57929 


.53832 


.49928 


8 


.72512 


.68946 


.65430 


9 


.83060 


. 8C257 


.77402 


10 


.90121 


. 88087 


.85947 


11 


.94542 


. 93169 


.91671 


12 


.97141 


. 96275 


.95292 


13 


.98579 


.98067 


.97460 


14 


.99328 


.99044 


.98691 



20 


25 


30 


.00047 


.00005 


.0 


.11016 


.05062 


.02283 


.45822 


.32060 


.22015 


.75491 


.63168 


.51860 


.91124 


.83568 


.75243 


.97403 


.93814 


.89000 


.99390 


.98034 


.95747 


.99885 


.99474 


.98570 


.99981 


.99881 


.99582 


.99995 


.99975 


.99892 

.99973 

.99991 


55 


60 


65 


.02914 


.01913 


.01252 


.16749 


. 13154 


.10281 


.38532 


.33157 


.28436 


. 59608 


.54125 


.48980 


.75745 


.71160 


.66627 


.86533 


.83164 


.79665 


.93062 


.90829 


.88387 


.96680 


.95332 


.93771 


.98524 


. 97779 


.96860 


.99389 


.99011 


.98512 


.99762 


.99587 


.99335 


.99911 


.99835 


.99718 


90 


95 


100 



00144 






02903 






12723 






28665 






46229 


.42740 


.39460 


61988 


.58639 


.55395 


74524 


.71647 


.68787 


83725 


.81443 


.79119 


90068 


. 88374 


.86605 


94204 


.93020 


.91749 


96764 


.95982 




98271 


.97781 





35 



• 0 

.01016 
.14921 
.42018 
.66858 
.83417 
.92640 
. 97099 
.98984 
.99683 
.99910 
.99974 



70 



.00817 

.08019 

.24317 

.44197 

.62203 

.76098 

.85782 

.92024 

.95777 

.97888 

.99001 

.99551 
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APPENDIX B. 









P< R 

N 


< K/N ) 






N = 


5 


10 


15 


20 


25 


30 


K 

I 


•27854 


.03814 


.00455 


.00051 


.00006 




2 


.83040 


.47338 


.23874 


.11374 


.05253 


.02377 


3 


.98464 


.83505 


.63716 


.46151 


.32425 


.22334 


A 


. 9997C 


.96661 


. 87390 


.75422 


.63247 


.52053 


5 


1.0 


.99582 


.96699 


.90962 


.83454 


.75201 


6 




.99970 


.99358 


.97307 


.93696 


.88891 


7 




.99998 


.99909 


.99354 


.97969 


.95663 


8 




.99999 


. 99989 


.99876 


.99448 


.98526 


9 




. 99999 


. 99997 


.99980 


.99873 


.99564 


10 




1.0 


.99998 


.99996 


.99974 


.99997 


N = 


35 


40 


45 


50 


55 


60 


K 

L. 


. 83336 


.77405 


. 71369 


.65414 


.59666 


. 54205 


7 


.92551 


. 88924 


. 34672 


.80260 


.75721 


.71152 


a 


. 97039 


.95010 


.92511 


. 89635 


.86476 


.83117 


9 


.98954 


.97988 


. 96658 


.94983 


.93006 


.90773 


1C 


.99671 


.99268 


.98638 


.97761 


.96639 


.95237 


11 


.99906 


. 99758 


.99492 


.99078 


.98499 


.97748 


12 


.99974 


.9^925 


. 99825 


.99648 


.99376 


.98993 


13 


.99990 


.99976 


.99942 


.99873 


.99756 


.99577 


14 


.99994 


.99990 


. 99979 


.99955 


.99908 


.99831 


N = 


65 


70 


75 


80 


35 


90 


K 

6 


.49076 


.44302 


.39892 


.35854 


.32142 


.28773 


7 


.66645 


. 62239 


.57979 


. 53893 


.49999 


.46306 


a 

W 


.79632 


.76078 


.72506 


.68953 


.65448 


.62016 


9 


.68337 


.85739 


. 93024 


.80229 


.77383 


.74513 


10 


.93725 


.91979 


.90077 


.88047 


.85911 


.83694 


11 


.96828 


.95714 


.94504 


.93130 


.91633 


.90031 


12 


.98492 


.97864 


. 97114 


.96244 


.95261 


.94171 


13 


.99324 


. 98986 


.98562 


.98046 


.97438 


.96739 


14 


.99713 


.99543 


.99316 


.99031 


.98677 


.98253 
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APPENDIX C. 


ASYMPTOTIC 


COMPARISON 


OF D 

N 


♦ C » 
N 


ANO R 

N 




N = 


1C 


15 


20 


25 


30 


35 


ALPHA 


.20 


A 


.322 


.266 


.231 


.210 


. 192 


. 181 




D 


.338 


.284 


.239 


.212 


. 196 


. 181 




C 


.292 


.246 


.205 


. 193 


. 178 


. 166 




P 


.309 


.245 


.205 


.193 


.178 


. 166 


.10 


A 


.368 


• 3C4 


.264 


.238 


.221 


. 206 




D 


.386 


.316 


.273 


.244 


.223 


. 206 




C 


• 348 


.284 


.246 


.226 


.205 


. 195 




P 


.342 


. 284 


.246 


.226 


.205 


. 195 


• C5 


A 


.409 


.338 


.294 


. 264 


.242 


. 230 




D 


.431 


. 352 


.304 


.272 


.248 


.230 




C 


.386 


. 321 


.281 


.251 


.230 


.217 




R 


.386 


. 321 


.281 


.251 


.230 


.217 


.025 


A 


.463 


. 378 


.327 


.293 


.268 


.248 




D 


.453 


. 364 


.324 


.292 


.268 


. 248 




C 


.425 


.351 


.300 


.275 


.254 


. 244 




P 


.428 


. 351 


.300 


.275 


.254 


.244 


.01 


A 


.516 


.423 


.365 


.326 


.298 


.276 




n 


.487 


.404 


. 352 


.321 


.293 


. 276 




c 


.478 


. 390 


.343 


.308 


.283 


. 258 




p 


.481 


. 390 


.343 


.308 


.283 


.258 





N= 


40 


45 


50 


60 


N>60 


ALPHA 

.20 


A 


.169 


.160 


.157 


• 138 






D 


.169 


. 160 


. 157 


.138 


1.07/lfF? 




C 


.156 


. 148 


. 140 


.134 




P 


.156 


. 148 


. 140 


. 134 


. 10 


A 


. 193 


.182 


.173 


. 158 






D 


.193 


.182 


.173 


. 158 


1.22 /"Yn 1 




C 


. 179 


.170 


.161 


.155 




P 


.179 


. 170 


. 161 


.155 


.05 


A 


.216 


.203 


.194 


. 176 






D 


.216 


.203 


.194 


. 176 


1.36/ir 




C 


. 200 


.190 


.180 


.173 




P 


. 200 


.190 


.180 


. 173 


.025 


A 


.231 


.218 


.207 








D 


.231 


.218 


.207 




1.43/1^ 




C 


. 222 


. 209 


.200 






P 


.222 


.209 


.200 




.01 


A 


. 258 


.243 










D 


.258 


. 243 






1.63/Vn 1 




C 


.249 


. 243 








P 


.249 


.243 








A = 


ASYMPTOTIC ESTI 


MAT I ON 







D = D 

N 

C = C 

N 

R = R 

N 
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COMPUTER OUTPUT 



SAMPLE PROGRAM 
THE D STATISTIC FOR N=10 



I A = 


1 


I B = 


1 


I C = 


0 


NNN= 


I A = 


2 


I B = 


2 


IC = 


0 


NNN = 


I A = 


3 


I B = 


3 


IC = 


0 


NNN= 


I A = 


4 


I B = 


4 


IC = 


0 


NNN= 


I A = 


5 


I B= 


5 


IC = 


0 


NNN= 


I A = 


6 


I B = 


6 


I c = 


0 


NNN= 


I A = 


7 


IP = 


7 


IC = 


0 


NNN= 


14 = 


8 


16 = 


8 


IC = 


0 


NNN= 


I A = 


9 


I B = 


9 


IC = 


0 


NNN= 



OUT OF RANGE IA= 10 IB= 
THE ABOVE IS AN EXAMPLE OF 



10 


PROB. = 


0.0003628 


10 


PROB. = 


0.2512797 


10 


PROB. = 


0.7294613 


10 


PROB. = 


0.9410061 


10 


PROB. = 


0.9922168 


10 


PROB. = 


0.9994260 


10 


PROB. = 


0.9999743 


10 


PROB. = 


0.9999936 


10 


PROB. = 


0.9999938 


10 


IC= 0 


NNN= 10 



ONE OF THE ERROR MESSAGES 
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COMPUTER PROGRAM 



1 FCRMAT(16X»* SAMPLE PROGRAM*, /) 

2 FORMAT ( 1 OX , * THE D STATISTIC FOR N=10*,/) 

3 FORMAT ( 1 X, ' THE ABOVE IS AN EXAMPLE OF*, 
1* ONE OF THE ERROR MESSAGES.*/) 

C 

C 

WR I TE ( 6 , 1 ) 

WPITE(6,2) 

IC=0 

NNN=10 

DC 10 I A=1 ,10 

10 CALL DURBINC I A , I A , IC , NNN ) 

WP I TE ( 6 , 3 ) 

STOP 

END 



SUBROUTINE DURB I N ( I A , I B , I C , NNN ) 

THIS SUBROUTINE COMPUTES THE PROBABILITY THAT A 
RANDOM SAMPLE OF SIZE NNN DRAWN FROM A UNIFORMS, 1) 
RANDOM VARIABLE AND PLOTTED AS SHOWN IN SECTION 
FOUR OF THIS THESIS WILL STAY WITHIN THE REGION, R, 
BETWEEN THE LINES, 

NNN*Y=NNN*IA+(NNN+IC)*X 
NNN*Y= -NNN*IB+( NNN+IC)*X 



THE RESTRICTIONS ARE 

IA, IB, ANC IC MUST BE INTEGER VARIABLES 
IC.GE.O 

IA. GT.0 

IB. GT.O 
IA+IC.GT.O 
IB-IC.GT.O 

3 FORMAT ( /// ,2X, ' OUT OF RANGE IA=*,I3,» IB=',« IC= *, 
113,* NNN= *,I3,//) 

4 FORMATt//,* IP.GT.NNN IA=*,I3,* IB=*,I3, 

1* IC = * » 13, • IP= * , 1 3 , ' NNN= • , I 3 ) 

5 FORMATt* I A= * , I 3 , ' I B= * * I 3 , • IC = *,I3,» NNN = *,I3, 

1 • PROB. = ' , F9 • 8 ) 



DIMENSION H( 1600) , HSTORE ( 1600 ) ,SAVE1( 1600) ,SAVE2( 1600) 
DATA H/16C0*0.0/ 

DATA HSTORE/ 1600*0.0/ 

DATA SAVE1 /1600*0.0/ 

DATA SAVE2/1600*0.0/ 

DATA JPH/O/ 

IF(IA.LE.O) GO TO 10000 
IF(IB.LE.O) GO TO 10000 
I F ( ( I A+ IC ) • LE. 0 ) GO TC 10000 
IF( C IB-IO.LE.C) GO TO 10000 
IF( JPH.GT.O) GO TO 100 
CALL DURMAT(HSTORE ,40 ) 

JPH=100 
100 IP=IA*IB+1 

I F ( IP.GT.4C) GO TO 10000 
IJ=IB*IP+IB-IC+1 

I F ( ( (IB+IC).GE.NNN).ANO. ( (IA+IO.GE.NNN) ) GO TO 11000 
NS I ZE=IP*IP 
JGOT =0 

DC 500 1=1 , IP 
JGE T= ( I- 1 ) *40 
DO 50C J*1,IP 
JGET=JGET+ 1 
JGOT=JGOT+ 1 
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500 H( JG0T)=HST0RE( JGET) 

DO 600 J=1,IP 
M=J*IP 

600 H(M)=0.0 

DO 700 J=1,IP 
M= ( J-l ) *1 P+J 
700 SAVE2 ( M ) =1 • 0 
I C K = NNN+ 1 C 

FACT=REDUCE(NNN, IC,1 ) 

750 CONTINUE 

DO 8CC J=1 »NSI ZE 
800 SAVE2(J)=SAVE2( J ) *F ACT 

CALL GMPRD( SAVE2 , H , SA VE 1 , IP, IP, IP) 
I P ATH= 1 
I CK=ICK— 1 

I F ( ICK.LE.O) GC TO 2000 
FACT=PEDUCE(NNN,IC,2) 

DO 900 J=1,NSIZE 
900 SAVE1( J)=SAVEI ( J) *FACT 

CALL GMPRD( SAVEl ,H,SAVE2, IP, IP, IP) 
I P ATH = 2 
I CK= IC K— 1 

I F ( ICK.LE.O) GO TO 2000 
FACT-R EDUCE ( NNN * I C » 2 ) 

GO TO 750 

2000 GO TO (2010,2100), IPATH 

2010 WRITE (6,5) I A, IB, IC,NNN,SAVE1 ( I J ) 

GO TO 2500 

2 IOC WRITE (6,5) I A , I B , IC ,NNN , SAVE 2 ( I J ) 
2500 DC 2600 I=1,NSIZE 
2600 SAVE2 ( I )=G. 0 
RETURN 

10000 CONTINUE 

WPITE(6»3) IA,IB,IC,NNN 
RETURN 

11000 WRITE(6,4) IA, I B, IC, I P,NNN 
RETURN 
END 



FUNCTION REDUCE ( NNN, I C,IENTER) 

THIS SUBROUTINE COMPUTES THE TRANSITION MATRIX 
H OF SECTION FOUR OF THIS THESIS FOR INTEGER VALUES 
OF IA, IB, AND IC. 



GO TO Cl. ICO ) , I ENTER 
1 CONTINUE 
Z=NNN+IC 
W=NNN 
ICK=NNN-1 
Y=0. 0 
I PATH=2 
SAVE=W/Z 
1 « 

REDUCE=SAVE 

RETURN 

100 I F ( ICK.LE.C) GO TO 10C0 
GO T0( 200 , 300 ) , IPATH 
200 I PATH=2 
I CK= ICK-1 
SAVE=W/Z 
W=W-1.0 
REDUCE=SAVE 
RETURN 
30C I PATH=1 
Y=Y+ 1.0 
I C K = I C K- 1 
REDUCE=Y/Z 
RETURN 

1000 REDUCE=1.0 
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RFTURN 

END 



SUBROUTINE GMPRD ( A , B , R , N , M, U 

C THIS SUBROUTINE COMPUTES THE MATRIX PROOUCT OF 

C A*B AND STORES THE RESULTS IN R. 

DIMENSION A(1)«B(1),R(1) 

IP=0 
I K=-M 

DO 10 K=1,L 

IK=IK4M 

DO 10 J=1,N 
IR=I R+ 1 
JI=J-N 
I B = I K 
R ( IR )=0.0 
DO 10 1=1, M 
J I = J I + N 
I B=I P+ 1 

10 R ( I R I = P ( I R ) 4- A ( J I >*B( IB) 

RETURN 

END 



SUBROUTINE DURM AT ( CUR , JCUR ) 

THIS SUBROUTINE PROVIDES THE SCALE FACTOR CF 
P(IAylBflC) OF SECTION FOUR IN THE FOLLOWING ORDER, 
N/N,1/N,(N-1 )/N,2/N,.,. 

DIMENSION DUR(l) 

NNN=JDUR*JCUR 
DO 10 1=1, NNN 
10 DUP ( I J =0 • C 
I J= JDUR 

DO 200 1=2, JDUR 
I J=I J+I-l 
DUR ( I J > =1 . 0 
I J=I J+l 
DUR ( I J 1 =1 . 0 

I F ( I J.GE.NNN) GO TO 3C0 
NN= I +1 
Z = 2.0 

DC 100 I K= NN , J DUR 
I J=I J+l 

DUR ( I J ) =DUR ( IJ-1)*( 1.0/Z) 

100 Z=Z+1.C 
200 CONTINUE 
300 RETURN 
END 
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